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ABSTRACT 

We establish new dynamical constraints on the mass and abundance of compact objects in the halo 
of dwarf spheroidal galaxies. In order to preserve kinematically cold the second peak of the Ursa Minor 
dwarf spheroidal (UMi dSph) against gravitational scattering, we place upper limits on the density of 
compact objects as a function of their assumed mass. The mass of the dark matter constituents cannot 
be larger than 10'^ Mq at a halo density in UMi's core of 0.35 Mq pc~'^. This constraint rules out 
a scenario in which dark halo cores are formed by two-body relaxation processes. Our bounds on the 
fraction of dark matter in compact objects with masses > 3 x 10^ M0 improve those based on dynamical 
arguments in the Galactic halo. In particular, objects with masses ~ 10^ Mq can comprise no more than 
a halo mass fraction ~ 0.01. Better determinations of the velocity dispersion of old overdense regions 
in dSphs may result in more stringent constraints on the mass of halo objects. For illustration, if the 
preliminary value of 0.5 km/s for the secondary peak of UMi is confirmed, compact objects with masses 
above ~ 100 Mq could be excluded from comprising all its dark matter halo. 

Subject headings: dark matter — galaxies: halos — galaxies: individual (Ursa Minor dwarf spheroidal) 
— galaxies: kinematics and dynamics — galaxies: structure 



1. INTRODUCTION 

The composition of dark halos around galaxies is a difh- 
cult problem. Many of the baryons in the universe are dark 
and at least some of these dark baryons could be in galac- 
tic halos in the form of very massive objects (VMOs), with 
masses above 100 Mq. Astrophysically motivated candi- 
dates include massive compact halo objects (MACHOs) 
and black holes either of intermediate mass (IMBHs; 10^ '^ 
to 10^ Mq) or massive (> 10^ Mq). IMBHs are an in- 
tringuing possibility as they could contribute, in principle, 
to all the baryonic dark matter and may be the engines 
behind ultraluminous X-ray sources recently discovered in 
nearby galaxies. 

A successful model in which VMOs are the dom- 
inant component of dark mat ter halos could r esolve 
some long-standing problems (iLacev fc Ostrikeil Il985t 
iTremaine fc Ostriked[l99i iJin et al.ll2005[r If the halos 
of dSphs are comprised by black holes of masses between 
^ 10^ and 10^ Mq, they evolve towards a shallower inner 
profile in less than a Hubble time, providing an explana- 
tion for the origin of dark matter cores in dwarf galaxies, 
and the orbits of globular cluste rs (GCs) do no t shrink to 
the center by dynamical friction (|Jin et al.ll2005r ). Very few 
observational limits on VMOs in dSph halos have been de- 
rived so far. This Letter is aimed at constraining the mass 
and abundance of VMOs in the halos of dSphs by the dis- 
ruptive effects they would have on GCs and cold long-lived 
substructures. 



2. CONSTRAINTS FROM THE SURVIVAL OF FORNAX GCS 

Fornax is a dark dominated dSph galaxy with an unusu- 
ally high GC frequency for its dynamical mass. In order 
to place dynamical constraints on the mass of VMOs by 
requiring that not too many GCs are disrupted, the den- 
sity of VMOs ph with mass Mh along the orbits of the 



GCs should be known. Since the three-dimensional dis- 
tances of the GCs to the center of Fornax are unknown, 
the density of VMOs at a mean distance of ~ 1 kpc will 
be adopted. By / we will denote the halo mass fraction 



in VMOs, i.e. / = ph/pdm, where pd 



0.02-0.05 Mq pc"^ a t 1 kpc (jWalker et al.ll2006, 



is the halo densit 



TlClessen fc BurkertI (|1996[ ) give the "survival diagram" 
of GCs, considering different encounter histories for ph = 
0.026 Mq pc"^ and a velocity dispersion of halo particles 
(Tfi = 120 km s~^. The survival diagram establishes the 
range of mass and concentration such that GCs with cen- 
tral densities between 10^ Mq pc~^ and 10^ Mq pc~^ 
(or, equivalently, core radii between 0.5 and 2 pc) have a 
probability of less than 1% to survive after 10 Gyr. Since 
the dissolution timescale for a certain GC (and a given 
Mil) scales as oc ah, and cr/i « 20 km/s in Fornax, it 
will be a factor of (5-12)/ less in Fornax as compared t o 
the Galactic case considered bv lKlessen fc BurkertJ(ll99fih . 
Therefore, the survival diagram for Fornax's GCs can be 
derived at once (see Fig. [T]). GCs in the region above the 
thick line have a probability of less than 1% to survive in 
a dark halo with phMh = lO'^ Mq pc~^. For such a puMh 
value and if the present parameters of the GCs are repre- 
sentative of their parameters at the time they formed, we 
would be observing the lucky survivors of an initial popu- 
lation of > 500 GCs of 2 X 10^ Mq, which turns out to 
be very unrealistic for a galaxy with a V-band luminosity 
of 1.5 X 10^ Lq. 

It is possible to estimate the probability that we 
are observing only the survivors of a larger original 
population that are in th e process of quick disruption 
(jTremaine fc Ostrikei1l999f) . The distribution of the disso- 
lution age of GCs is expected to follow a scale-free power 
law F = C(tdis + tH)~'^, where idis is the characteristic 
dissoluti on timescale and tn is th e age of the cluster pop- 
ulation ([GnedinXioitrikeB [l993) . Fig. [1] shows that at 
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least 4 GCs are above the thick line; this indicates that 
their present disruption timescales are < 0.22tH- If the 
exponent of distribution of lifetimes i s q ^ 2, as derived 
for Galactic GCs (,Gnedin fc Ostriker 1997f ). the probabil- 
ity to have 4 out of 5 GCs with lifetimes less than Q.22tH 
is ^ 1%, whereas the probability that the lifetimes of all 
the GCs are less then O.btn is 0.4%. Hence the probabil- 
ity that the whole dark halo is comprised of objects with 
masses > 5 x 10^(^,^/0.02 Mgpc-^)-! Mq is less than 1%. 
If only a fraction / of the dark mass is in compact objects 
of mass Mh, then / < 5 x 10"^ Mq/M^. 

3. PERSISTENCE OF DYNAMICALLY-COLD 
SUBPOPULATIONS 

Localized regions with enhanced stellar density and, 
where data permit, extremel y cold kinematics have been 
detect ed in some d Sphs f e.g.. [Olszews ki fc Aronson' 1985': 
Klevna et all |2003| . hereafter K03; [Coleman et al, ,2004; 
Walker et al . 2006b). In particular, UMi dSph has re- 
ceived the most attention. Collecting the velocity of stars 
in 6' radius aperture, K03 found that a two-Gaussian pop- 
ulations, one representing the underlying 8.8 km/s Gaus- 
sian and the other with velocity dispersion CTs — 0.5 km/s, 
representing a subpopulation of fraction 0.7, is > 3 x 10^ 
times more likely than the default 8.8 km/s. The best-fit 
as is ill-determined as it is much smaller than the median 
velocity errors (5 km/s). Nevertheless, even with these 
fiducial errors, we can be certain that the velocity dis- 
persion is < 2.5 km/s at ~ 95% confidence level. The 
stars which form the secondary density peak are not dis- 
tinguished in colour and magnitude from the remainder of 
the UMi population (jKlevna et al.ll998D . In fact, UMi star 
formation history indicates that its stars ha ve been formed 
in a single burst earlier than 10 Gyr ago (jCarrera et al.l 
12002'). 

Although UMi has long been suspected of experienc- 
ing ongoing tidal disruption, regions with enhanced vol- 
ume density and cold kinematics cannot be the result 
of tidal interactions because the coarse-grained phase- 
density, p/cr^, in coUisionless sytems must be constant 
or even decrease, thus implying that overdensity regions 
should appear dynamically hotter. This suggests that the 
clump is long-lived. An alternative explanation is that 
the density peak is a projection effect and that what we 
are seeing is a cold, low-density tidal tail. However, nu- 
merical e xperiments have s hown that this scenario is very 
unlikely (|Read et al.ll2006[ ). The most plausible interpre- 
tation is that the clump is a disrupted stellar cluster, now 
surviving in phase-space because the underlying gravita- 
tional potential is harmonic (K03). Within this potential, 
gravitational encounters with the hypothetical VMOs will 
dominate the orbital diffusion of the stars once they be- 
come unbound from the progenitor cluster. The integrity 
of the cold clump may impose useful upper limits on the 
mass of VMOs. The fact that the subpopulation is orbiting 
within the dark matter core of UMi will greatly simplify 
its dynamical description. 

Clump's stars will undergo a random walk in momen- 
tum space by the collisions with the population of VMOs. 
Here we are interested in the velocity change induced in a 
star relative to the clump center of mass. The mean-square 
velocity change of a star in an encounter with a VMO of 



mass Mh, and impact parameter b > 5ri/2, with ri/2 the 
clump's median radius, in the impulsive approximation is: 



- 2 f 2GAh Y^ 



(1) 



where V is the maximum relative velocity between the 
clump and the perturber and is the rnean-s quare po- 
sition of the stars in the clump (jSpitzed Il958f ). In the 
opposite case of a head-on collision {b = 0), the mean 
change is comparable to that predicted by the tidal ap- 
proximation when b ~ 1.4ri/2. The usual way to proceed 
is to integrate Av"^ given in Eq. ([T]) for impact parameters 
b from 1.4ri/2 to infinity and correct for the encounters 
in w hich the tidal approximation fails by a facto r g ~ 3 
fe.g.. lBinnev fc Tremainlll987t iGieles et al.ll2006f ). Doing 
so, and for a distribution of clumps and VMOs with a rel- 
ative one-dimensional velocity dispersion firci, we obtain: 



(2) 



For UMi, the persistence of the clump for a large frac- 
tion of a Hubble time indicates a core of the dark halo 
of at least 2-3 times the size of the orbit of the clump, 
which is > 150 pc. In terms of the stellar core radius 
(~ 200 pc), this makes a halo core 1.5-2 times the stellar 
core and, consequently, the velocity dispersion of the halo 
particles in the core is, at least ~ 1.5-2 times the stellar 
velocity dispersion, corresponding to Uh ^ 15-20 km/s. 
The impulsive approximation is valid for bu < V , where 
oj — (Ts/f 1/2 ^-nd Us is the internal one-dimensional velocity 
dispersion of the subpopulation. For the encounters with 
b < 5r]^/2j responsible for most of the velocity impulse, 
this condition is well satisfied for ah ^ has- Therefore, 
for as ''^ 1 km/s, this requirement is fulfilled within the 
isothermal dark core of UMi. 

Since stars in the clump are unbound, the self-gravity 
of the clump in a first approximation can be ignored con- 
sidering only orbit diffusion in the large-scale harmonic 
potential of the parent galaxy. In a one-dimensional har- 
monic potential, a velocity impulse Av'^ produces a change 
in the velocity dispersion Aa^ = A (w'^) — Av^ /2, where 
the brakets (...) refer to the mean value after averaging 
over one orbit. Combining this relation with Eq. ([2]), we 
find the change of as in a time At: 



Aai 



S^gG^PhA'hr^At 



27ahrl 



(3) 



/2 



where cTrci ~ ah is assumed, since the population of clumps 
is expected to have a velocity dispersion similar to the stel- 
lar background, ~ 9 km/s in UMi. The ratio i] = 
depends on the model: for a Plummer cluster rj — A, 
whereas 77 = 1.5 for both a King pr ofile with a dimen sion- 
less potential depth of Wq — 9 (e.g.. lGieles et al.ll2006D and 
a Gaussian density distribution. In order to take a conser- 
vative value and to facilitate comparison with photometric 
and theoretical analysis that assume Gaussian models, we 
will adopt ?7 = 1.5 constant in time. 

By i2.5 we will indicate the time required for a very cold 
group of unbound stars as — 0, to acquire a velocity dis- 
persion of 2.5 km/s. From Eq. ^ with g = 3 we then 



3 



obtain: 



'-2.5 



5Gyr 



Ph 



O.lM0/pc3 



Mh 



5 X 103 Mf; 



20 km/s 

Since many of the recent dynamical models suggest 
mean dark matter densities within the stellar core ra- 
dius of UMi > 0.1 Mq pc"3j correspondin g to a cen- 
tral mass-to-Ught ratio of > 15 M^/Lp) (iLake 
Prvor fc Kormendyl Il99d; llrwin &: Hatzidimitriou 



Mated 11998: Wilkinson ct al. 200"^ 



199C; 



1995 



0.1 Mq pC" 



will be accepted as a conservative reference value. 

The mean-square radius of the subpopulation increases 
in time with as, according to the relation ~ cr^/il^, 
where fl is the orbital frequency in the constant-density 
core. If the core is dominated by the dark matter compo- 
nent pdm/^'^ = {4:TtG/3)-\ then 

^ = /^i. (5) 

VTT ah 



If as = 2.5 km/s, the mean radius vr^ ~ ag/^ ~ 60 pc. 
This corresponds to a 1 cr radius of 25.5 pc for a Gaussian 
density profile, implying an angular size of almost 2' at 
the distance of UMi (~ 66 kpc) . This value is comparable 
to but slightly larger than the observed 1 a radius of the 
secondary peak ~ 1.6'. As a consequence as < 2.5 km/s; 
otherwise, the stellar subpopulation would appear more 
extended and diffuse than it is observed. This upper value 
is in agreement with the observations of K03. 

Our estimates for the size and velocity dispersion of the 
subpopulation are independent of the eccentricity of the 
orbit because there is little variation of the macroscopic 
properties of the halo, ph, ah or Q, within the core where 
the clump is orbiting. Nevertheless, the dark halo could 
have suffered significant evolution due to two-body pro- 
cesses and tidal stirring. For Mh < 10^ Mq relaxation 
processes induce an insig nificant change in the internal 
properties of the halo fe.g.. lJin et al.|[2005[ ). Tidal heating 
can lead to a reduction of both the density of dark mat- 
ter particles and its vel ocity dispersion fe.g.. iMaver et all 
[2OOII: iRead et al.l [20061) . Since the phase-space density, 
Ph/cr\, for coUisionless systems is nearly constant or de- 
creases with time, the rate of energy gained by the clump 
due to encounters with VMOs should have been more in- 
tense in the past. The inclusion of evolution of halo prop- 
erties by tidal effects would lead to a stringent upper mass 
limit. 

If the stellar progenitor cluster became unbound imme- 
diate ly after formatio n when supernovae expel the gas con- 
tent (|Goodwinlll997t K03), ^2.5 should be greater than the 
age of the cluster tn ~ 10 Gyr. This requirement com- 
bined with Eq. (jU implies: 

A4<2.5x103Mo('— -^V^) '(-^\. (6) 



,0.1Mq/pcV V20km/s, 
If initially the stellar cluster is gravitationally bound, 
the increase in the internal energy gained by encounters 
with VMOs will eventually exceed its binding energy af- 
ter a time ibe- Let us estimate the, the time at which 
the cluster becomes unbound. K03 infer a total mass 
of the cluster, M^, of 3 x lO" Mq. If this cluster fol- 
lowed t he recently ob served relation between radius and 
mass of iLarseiil ()2004[ ). ri/2 « 4 pc. Adopting the refer- 
ence values of the dark matter halo (pdm = 0.1 M© pc"^ 



and ah = 20 km/s) and rescaling the survival diagram of 
iKlessen fc BurkertI (jl996[ ) (their figure 11) for the param- 
eters of UMi, we infer that more than 95% of the clusters 
with mass 3 x 10^ M© will become unbound after the ~ 3 
Gyr if Mh > 3.5 x 10^/-^ Mq. Therefore, in order to have 
a dynamically cold subpopulation with as < 2.5 km/s, as 
that observed in UMi at the present time tn ~ 10 Gyr, 
we need i2.5 > iff ~ ^fce, which implies the following upper 
limit for Mh 

\ -1 / 

-,3^;r / Ph \ I <^h 



Mh < 3.5 X 10^ Me 



(7) 



^O.IMq/pcV V 20 km/s. 
Other corrosive effects such as mass loss by stellar evolu- 
tion or tidal heating may also accelerate the dissolution of 
the cluster. Therefore, our estimates for the and, hence, 
for Mh, are upper limits. 

The survival probability after collisions with VMOs in- 
creases for progenitors that are more compact. For in- 
stance, if the progenitor were a supercluster with a core 
radius Vc ~ 0.5 pc and central density ^ 3 x 10*^ Mq pc~^, 
the probability of its remaining gravitationally bound after 
6 Gyr is ~ 25%, for Mh = 6.5 x 10^ f'^ Mq. Hence, there 
may be a non-negligible probability that such a superclus- 
ter has survived bound for 6 Gyr and that during the sub- 
sequent 4 Gyr it is dynamically heated by 6.5 x 10^/"^ 
M© VMOs to reach as = 2.5 km/s at the present time. 
Unfortunately, the evaporation time for this supercluster, 
setting the scale for dynamical dissolution by internal pro- 
cesses, is very short. In fact, for such a stellar cluster, 
the evaporation time is ^ 20^^.^, with tr h the half-mass 
relaxation time (jGnedin fc Ostrikei1ll997l ). The resulting 
evaporation time is < 1 Gyr for an average stellar mass 
> 1 Mq and, hence, internal processes would have pro- 
duced a fast desintegration of such a cluster. We conclude 
that the upper limit given in Eq. ([7|) is realistic and robust. 

Our approximations break down when the halo only con- 
tains a few VMOs; at least 5 objects within a radius of 600 
pc are required, implying that our analysis is restricted to 
masses Mh < 2 x 10''/ Mq. In Fig. ^ the observational 
limits on VMOs over a wide range of masses and dark 
matter fractions are shown. 

4. DISCUSSION AND CONCLUSIONS 

The analysis of the survival of Fornax's GCs rules out 
the mass range that would be interesting for explaining 
the origin of dark matter cores in dwarf galax ies (K03; 
iGoerdt et aT]|2006l : Isillchez-Salcedo et al.ll2006[ ). because 
the relaxation timescale for VMOs of mass < 5 x 10^ Mq 
exceeds the Hubble time. Moreover, it was found that 
the integrity of cold small-scale clustering seen in some 
dSphs imposes more stringent constraints on the mass of 
VMOs. A source of uncertainty is the mean density of 
dark matter within the core of dSphs. In the particular 
case of UMi and according to the scaling relations com- 
piled in iKormendv fc FreemanI (j2004,), the corresponding 
central dark matter density is 0.35 Mq pc"^. A sHghtly 
larg er value has been d erived from its internal dynam- 
ics (jWilkinson et al.l [2006^) . At a density pdm — 0.35 
Mq pc-3, Eq. implies that Mh < 1000/"! Mq. We 
strongly encourage better determinations of the velocity 
dispersion of cold density aggregates (boimd or unbound) 
in dSphs. For instance, if the preliminary quoted value of 
0.5 km/s for the secondary peak of UMi were confirmed. 
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our upper limit for Mh would be immediately reduced by 
a factor of 25, implying a very tight bound Mh < (40- 
120)/"^ Mq, depending on the adopted dark matter den- 
sity (0.3 Mq pc~3 to 0.1 Mq pc-3). 

In a unifi e d sch eme such as the 'stirring scenario' by 
iMaver et al.l (|200lh . the composition of the dark halos of 
low-s urface brightness a nd dSph galaxies should be the 
same. iRix fc Lakd (L1993) found an upper limit of lO'' Mq 
by examini ng the dynamical heating o f the stellar disk 
of GR iTremaine fc Ostriked (Il999l) warn a bout the 
weakness of the argument of IRix &: Lakd (fl993 h because 
the rotation curve of GR 8 decays in a Keplerian fashion 
suggesting that GR 8 does not host a halo of dark matter, 
as expected if the dark halo has evaporated by two-body 
collisions. However, there exist dw arf galaxies with flat ro - 
tation curves. From the sample of iHidal go-Game^ (j2004l ). 
we have selected dwarf galaxies with high inclination an- 
gles and estimated the maximum permitted value of 
consistent with the thickness of the old stellar disk. Per- 
haps one of the most pristine cases is the edge-on galaxy 
UGCA 442. This galaxy shows the typical flat rotation 
curve, implying the existence of a halo of dark matter with 
a central density of 0.07 Mq pc~^ and a velocity dispersion 
ah = Vc/V2 ~ 35 km/s (|C6te et al.ll2000[ ). The vertical 



stellar velocity dispersion has been derived using vertical 
hydrostatic equilibrium cr^ = ttGE/i, with E the total cen- 
tral surface density of the disk 65 Mq pc~^ and h the 
scale height of the old disk derived photometrically to be 
350 pc. This yields az = 18 km/s. By requiring that 
Aaz < (Jz for At > 1.5 Gyr, which is the characteristic 
age of the old stellar population, the disk-heating argu- 
ment establishes Mh < 8 x 10* f^^ Mq (jLacev fc Ostrikeij 
,1985 ). 

Fig- m contains the most relevant dynamical constraints 
on VMOs. In the conservative case as = 2.5 km/s, the 
most stringent bound for masses between 10^ Mq and 10^ 
Mq comes from studies of the distributio n of separations o f 
wide stellar binaries in the Galactic halo (jYoo et al.ll2004l ). 
Still, masses in that range are permitted if halo objects are 
slightly extended (sizes > 0.05 pc). Other sources of un- 
certainty in this approach are the orbital distribution of 
the binaries and the duratio n of perturbations they are 
subjected to (jJin et al.ll2005[) . 
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Fig. 1. — Survival diagram for GCs in Fornax's halo for PhAIh = 10^ Mg pc ^. If GCs had core radii in the range from 0.5 to 2 pc, their 
probability of survival after 10 Gyr would be less than 1% in the region left and above the thick li ne. The core radii of Fornax 3, 4 and 5 lie 
within the mentioned range, but Fornax 1 and 2 present core radii of 10 and 5.6 pc, respectively IIMackev &: Gilmorc 2003) . Therefore, the 
probability of survival for the latter GCs will be much less than 1%. 
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VMO mass (Mq) 

Fig. 2. — Observational constraints on VMOs from MACHO microlcnsing experiments, the distribution of wide binaries in the thoroughly 
validated Bassel mass model of the Milky Way (Bissantz. Enjjlmaicr & Gerhard 2003), the heating of the Galactic disk, the heating of the 
stellar disk of UGCA 442, and the survival of UMi's dynamical fossil in UMi with = 0.35 Mq pc~'^. 



